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In this paper we show that the small world and weak ties phenomena can spontaneously 
emerge in a social network of interacting agents. This dynamics is simulated in the frame- 
work of a simplified model of opinion diffusion in an evolving social network where agents 
are made to interact, possibly update their beliefs and modify the social relationships 
according to the opinion exchange. 
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1. Introduction 

Modeling social phenomena represents a major challenge that has in recent years 
attracted a growing interest. Insight into the problem can be gained by resorting, 
among others, to the so called Agent Based Models, an approach that is well suited 
to bridge the gap between hypotheses concerning the microscopic behavior of indi- 
vidual agents and the emergence of collective phenomena in a population composed 
of many interacting heterogeneous entities. 

Constructing sound models deputed to return a reasonable approximation of 
the scrutinized dynamics is a delicate operation, given the degree of arbitrariness in 
assigning the rules that govern mutual interactions. In the vast majority of cases, 
data are scarce and do not sufficiently constrain the model, hence the provided 
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answers can be questionable. Despite this intrinsic limitation, it is however impor- 
tant to inspect the emerging dynamical properties of abstract models, formulated 
so to incorporate the main distinctive traits of a social interaction scheme. In this 
paper we aim at discussing one of such models, by combining analytical and nu- 
merical techniques. In particular, we will focus on characterizing the evolution of 
the underlying social network in terms of dynamical indicators. 

It is nowadays well accepted that several social groups display two main features: 
the small world property [18j and the presence of weak ties [12) . The first property 
implies that the network exhibits clear tendency to organize in large, densely con- 
nected, clusters. As an example, the probability that two friends of mine are also, 
and independently, friends to each other is large. Moreover, the shortest path be- 
tween two generic individuals is small as compared to the analogous distance com- 
puted for a random network made of the same number of individuals and inter-links 
connections. This observation signals the existence of short cuts in the social tissue. 
The second property is related to the cohesion of the group which is mediated by 
small groups of well tied elements, that are conversely weakly connected to other 
groups. The skeleton of a social community is hence a hierarchy of subgroups. 

A natural question arise on the ubiquity of the aforementioned peculariar as- 
pects, distinctive traits of a real social networks: how can they eventually emerge, 
starting from an finite group of randomly connected actors? We here provide an an- 
swer to this question in the framework of a minimalistic opinion dynamics model, 
which exploit an underlying substrate where opinions can flow. More specifically, 
the network that defines the topological structure is imagined to evolve, coupled 
to the opinions and following a specific set of rules: once two agents reach a com- 
promise and share a common opinion, they also increase their mutual degree of 
acquaintance, so strengthing the reciprocal link. In this respect, the model that 
we are shortly going to introduce hypothesize a co-evolution of opinions and social 
structure, in the spirit of a genuine adaptive network [131 119] . 

Working within this framework, we will show that an initially generated random 
group, with respect to both opinion and social ties, can evolve towards a final state 
where small worlds and weak ties effects are indeed present. The results of this 
paper constitute the natural follow up of a series of papers [3j El IS] , where the 
time evolution of the opinions and affinity, together with the fragmentation vs. 
polarization phenomena, have been discussed. 

Different continuous opinion dynamics models have been presented in literature, 
see for instance [TIB [IT], dealing with the general consensus problem. The aim is to 
shed light onto the assumptions that can eventually yield to fixation, a final mono- 
clustered configuration where all agents share the same belief, starting from an 
initial condition where the inspected population is instead fragmented into several 
groups. In doing so, and in most cases, a fixed network of interactions is a priori 
imposed [2], and the polarization dynamics studied under the constraint of the 
imposed topology. At variance, and as previously remarked, we will instead allow 
the underlying network to dynamically adjust in time, so modifying its initially 
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imposed characteristics. Let us start by revisiting the main ingredients of the model. 
A more detailed account can be found in [3]. 

Consider a closed group of N agents, each one possessing its own opinion on 
a given subject. We here represent the opinion of element i as a continuous real 
variable Oi G [0,1]. Each agent is also characterized by its affinity score with respect 
to the remaining N — 1 agents, namely a vector otij, whose entries are real number 
defined in the interval [0, 1]: the larger the value of the affinity Qy, the more reliable 
the relation of i with the end node j. 

Both opinion and affinity evolve in time because of binary encounters between 
agents. It is likely that more interactions can potentially occur among individuals 
that are more affine, as defined by the preceding indicator, or that share a close 
opinion on a debated subject. Mathematically, these requirements can be accom- 
modated for by favoring the encounters between agents that minimizes the social 
metric D\ 3 = |AO^|(l - a^) + A/}(0,cr), where AO^ = 0\ - O) is the opinions' 
difference of agents i and j at time t, and the last term is a stochastic contribution, 
normally distributed with zero mean and variance a. For a more detailed analysis 
on the interpretation of a as a social temperature responsible of a increased mixing 
ability of the population, we refer to [5J[SJ[S]. 

Once two agents are selected for interaction they possibly update their opinions 
(if they are affine enough) and/or change their affinities (if they have close enough 
opinions) , following: 



Fi {x) = tanrHPUX_ a c)) + ^ ^ ^ = _ _ Acg) (2) 



two activating functions which formally reduce to step functions for large enough 
the values of the parameters /3i and /32, as it is the case in the numerical simulations 
reported below. 

Despite its simplicity the model exhibits an highly non linear dependence on the 
involved parameters, a c , AO c and a, with a phase transition between a polarized 
and fragmented dynamics [3]. 

A typical run for N = 100 agents is reported in the main panel of Fig. [TJ for a 
choice of the parameters which yields to a consensus state. The insets represent three 
successive time snapshots of the underlying social network: The N nodes are the 
individuals, while the links are assigned based on the associated values of the affinity. 
The figures respectively refer to a relatively early stage of the evolution t = 1000, 
to an intermediate time t = 5000 and to the convergence time T c = 10763. Time 
is here calculated as the number of iterations (not normalized with respect to N). 
The corresponding three networks can be characterized using standard topological 
indicators [TJ [5] (see Table [JJ), e.g. the mean degree < k >, the network clustering 
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coefficient C and the average shortest path < £ >. An explicit definition of those 
quantities will be given below. 

In the forthcoming discussion we will focus on the evolution of the network 
topology, limited to a choice of the parameters that yield to a final mono cluster. 



Tabic 1. Topological indicators of the so- 
cial networks presented in Fig. [T] The 
mean degree < k >, the network cluster- 
ing C and the average shortest path < I > 
are reported for the three time configura- 
tions depicted in the figure. 





<k> (t) 


C(t) 


<£>(*) 


t = 1000 


0.073 


0.120 


3.292 


t = 5000 


0.244 


0.337 


2.013 


t = T c 


0.772 


0.594 


1.228 




5000 10000 15000 

time 



Fig. 1. Opinions as function of time. The run refers to a c = 0.5, AO c = 0.5 and a = 0.01. The 
underlying network is displayed at different times, testifying on its natural tendency to evolve 
towards a single cluster of affine individuals. Initial opinions are uniformly distributed in the 
interval [0, 1], while a? are randomly assigned in [0, 1/2] with uniform distribution. 
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2. The social network 

The affinity matrix drives the interaction via the selection mechanism. It hence can 
be interpreted as the adjacency matrix of the underlying social network, i.e. the 
network of social ties that influences the exchange of opinions between acquain- 
tances, as mediated by the encounters. Because the affinity is a dynamical variable 
of the model, we are actually focusing on an adaptive social network [13l [19] : The 
network topology influences in turn the dynamics of opinions, this latter providing 
a feedback on the network itself and so modifying its topology. In other words, 
the evolution of the topology is inherent to the dynamics of the model because of 
the proposed self-consistent formulation and not imposed a priori as an additional, 
external ingredient, (as e.g. rewire and/or add/remove links according to a given 
probability [TH [TS] once the state variables have been updated) . 

Remark 1. (Weighted network) Let us observe that the affinity assumes pos- 
itive real values, hence we can consider a weighted social networks, where agents 
weigh the relationships. Alternatively, one can introduce a cut-off parameter, uf. 
agents i and j are socially linked if and only if the recorded relative affinity is large 
enough, meaning ay > a/ . Roughly, the agent chooses its closest friends among all 
his neighbors. 

The first approach avoids the introduction of non-smooth functions and it is 
suitable to carry on the analytical calculations. The latter results more straightfor- 
ward for numerical oriented applications. 

As anticipated, we are thus interested in analyzing the model, for a specific 
choice of the parameters, a c , a and AO c , yielding to consensus, and studying the 
evolution of the network topology, here analyzed via standard network indicators: 
the average value of weighted degree, the cluster coefficient and the averaged short- 
est path. These quantities will be quantified for (i) a fixed population, monitoring 
their time dependence; (ii) as a function of the population size, photographing the 
dynamics at convergence, namely when consensus has been reached. 

2.1. Time evolution of weighted degree 

The simplest and the most intensively studied one-vertex (i.e. local) characteristic 
is the node degree B the total number of its connections or its nearest neighbors. 
Because we are dealing with a weighted network we can also introduce the weighted 
node degree, also called node strength [4j, namely Si(t) = ■ a\j/{N — 1). Its mean 

a Let us observe that the affinity may not be symmetric and thus the inspected social network will 
be directed. One has thus to distinguish between In-degree, ki n , being the number of incoming 
edges of a vertex and Out-degree, k ou t, being the number of its outgoing edges. In the following 
we will be interested only in the outgoing degree, from here on simply referred to as to degree. 
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value averaged over the whole network reads: 

1 - 

< s > (*) = ]y 5>*W • (3) 
i=i 

Let us observe that the normalization factor N — 1 holds for a population of N 
agents, self-interaction being disregarded, < s > belongs hence to the interval [0, 1] 
and having eliminated the relic of the population size, one can properly compare 
quantities calculated for networks made of different number of agents. 

All these quantities evolve in time because of the dynamics of the opinions 
and/or affinities. Passing to continuous time and using the second relation of ([TJ, 
we obtain: 

d 1 A d 



< s >- 



dt N(N - 1) .4^ dt 

Let us observe that the evolution of affinity and opinion can be decoupled when 
AO c = 1. For AO c < 1, this is not formally true. However on can argue for an 
approximated strategy [5], by replacing the step function T 2 by its time average 
counterpart 72, where the dependence in AO^ has been silenced. In this way, we 
obtain form ([4| 

d N 

di <S>= N(N- 1) S ~ 4) = T2(< s > - < s 2 >) , (5) 

ij"=l 

where < s 2 >= J2ct 2 j/(N(N — 1)). Let us observe that 72 is of the order of 1/N 2 
times, a factor taking care of the asynchronous dynamics [9]. 

In [6] authors proved that ([5]) can be analytically solved once we provide the 
initial distribution of node strengths (see |Appendix A for a short discussion of the 



involved methods). Assuming s,(0) to be uniformly distributed in [0,1/2], we get 
the following solution (see Fig. [2]): 

<S>{t) = ^—l-J^W l ^{-^) ' (6) 

Using similar ideas we can prove [6] that the variance a 2 (t) =< s 2 > — < s > 2 
is analytically given by 



2e 272 * 4e 272 * 

a 2 (t) 



log ' 



(e-»*-l) 2 (e-»* + l) (e^*-l) 4 L" b V 2 )\ ' ^ 
The comparison between analytical and numerical profiles is enclosed in Fig. [5J 
where the evolution of < s > (t) is traced. Let us observe that here 72 will serve 
as a fitting parameter, when testing the adequacy of the proposed analytical curves 
versus direct simulations, instead of using its computed numerical value [S]. The 
qualitative correspondence is rather satisfying, so confirming the correctness of the 
analytical results reported above. 
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Fig. 2. Evolution of < s > (t). Dashed line (blue on-line) refers to numerical simulations with 
parameters ct c = 0.5, AO c = 0.5 and <r = 0.3. The full line (black on-line) is the analytical 
solution (O with a best fitted parameter 72 = 1.6 10 -4 . The dot denotes the convergence time in 
the opinion space to the consensus state, for the used parameters affinities did not yet converge. 
Let us observe in fact that affinities and opinions do converge on different time scale [9]. 



Assume T c to label the time needed for the consensus to be reached. Clearly, 
T c depends on the size of the simulated system 0. From the above relation ([6]), the 
average node strength at convergence as an implicit function of the population size 
N reads: 



< a > (T c (Af)) 



<y 2 (N)T c (N) 



2 e 72(W)T c (JV) 



e 72(N)T c (N) _ 1 ( e 72(A0T c (A0 _ 1): 



■log 



(8) 

where we emphasized the dependence of 72 and T c on N. However, as already 
observed l2 (N) = O (iV~ 2 ) and T C (N) = 0(N a ), with a € (1,2). Hence 
72(A0r c (A0 -> when N -> 00 and thus < s > (T C (N)) is predicted to be a 
decreasing function of the population size N, which converges to the asymptotic 
value 1/4, a value identical to the initial average node strength (see Fig. [3j, given 
the selected initial condition. In sociological terms this means that even when con- 
sensus is achieved the larger the group the smaller, on average, the number of local 



b In [3|[Zl it was shown that T c scales faster than linearly but slower than quadratically with respect 
to the population size N. 
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acquaintances. This is a second conclusion that one can reach on the basis of the 
above analytical developments. 



as - 

0.7- 

V 



'5 55 105 155 205 255 305 355 405 455 505 550 
N 



Fig. 3. Average node strength at convergence as a function of the population size. Parameters 
are AO c = 0.5, cr = 0.5 and four values of a c have been used : (0) etc = 0, (A) a c = 0.25, (□) 
a c = 0.5 and (Q) a c = 0.75. Vertical bars are standard deviations computed over 10 replicas of 
the numerical simulation using the same initial conditions. 



2.2. Small world 

Several social networks exhibit the remarkable property that one can reach an arbi- 
trary far member of the community, via a relatively small number of intermediate 
acquaintances. This holds true irrespectively of the size of the underlying network. 
Experiments |16) have been devised to quantify the "degree of separation" in real 
system, and such phenomenon is nowadays termed the "small world" effect, also 
referred to as the "six degree of separation" . 

On the other hand several, models have been proposed [El [17] to construct 
complex networks with the small world property. Mathematically, one requires that 
the average shortest path grows at most logarithmic with respect to the network size, 
while the network still displays a large clustering coefficient. Namely, the network 
has an average shortest path comparable to that of a random network, with the same 
number of nodes and links, while the clustering coefficient is instead significantly 
larger. 
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On the evolution of a social network 9 

In this section we present numerical results aimed at describing the time evo- 
lution of both the average shortest path and the clustering coefficient of the social 
network emerging from the model. As before, the parameters are set so to induce 
the convergence to a consensus state in the opinion space. 

We will be particularly interested in their asymptotic solutions, terming the 
associated values respectively < I > (T c ) and C(T C ) once the consensus state has 
been achieved. 

In Fig. |4] we report these quantities (normalized to the homologous values esti- 
mated for a random network with identical number of nodes and links) versus the 
system size. The (normalized) clustering coefficient is sensibly larger than one, this 
effect being more pronounced the smaller the value of a c . On the other hand the 
(normalized) average shortest path is always very close to 1. 

Based on the above we are hence brought to conclude that the social network 
emerging from the opinion exchanges, has the small world property. This is a re- 
markable feature because the social network evolves guided by the opinions, as it 
does in reality, and not result from an artificially imposed recipe. 
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Fig. 4. Normalized clustering coefficient (left panel) and normalized average mean path (right 
panel) as functions of the network size at the convergence time. Parameters are AO c = 0.5, 
a = 0.5 and four values of ct c have been used : (0) ct c = 0, (A) ct c = 0.25, (□) a c = 0.5 and (O) 
etc = 0.75. Vertical bars are standard deviations computed over 10 repetitions. 
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2.3. Weak ties 

Social networks are characterized by the presence of hierarchies of well tied small 
groups of acquaintances, that are possibly linked to other such groups via "weak 
ties" . According to Granovetter [12] these weak links are fundamental for the cohe- 
sion of the society, being at the basis of the social tissue, so motivating the statement 
"the strength of weak ties" . 

The smallest group in a social network is composed by three individuals sharing 
high mutual affinities, in term of network theory they form a clique [T], i.e. a maximal 
complete graph composed by three nodes. This can of course be generalized to larger 
maximal complete graphs, defining thus m-cliques. 

The degree of cliqueness of a social network is hence a measure of its cohe- 
sion/fragmentation: the presence of a large number of m-cliques together with very 
few, m'-cliques, for m' > m, means that the population is actually fragmented into 
small pieces, of size m not strongly interacting each other. 

We are interested in studying such phenomenon within the social network emerg- 
ing from the opinion dynamics model here considered, still operating in consensus 
regime. To this end we proceed as follows. We introduce a cut-off parameter a/ 
used to binarize the affinity matrix, which hence transforms into a an adjacency 
matrix a. More precisely, agents i and j will be connected, i.e. Oy = 1, if and only 
if otij > a.f. Once the adjacency matrix is being constructed, we compute the num- 
ber of m-cliques in the network. Let us observe that this last step is highly time 
consuming, being the clique problem NP-complete. We thus restrict our analysis to 
the cases m e {3, 4, 5}. 

For small values of a/ the network is almost complete, while for large ones it can 
in principle fragment into a vast number of finite small groups of agents. As reported 
in the inset of right panel of Fig. [5j for a/ ~ 1 only 3-cliques are present. Their 
number rapidly increases as a/ is lowered. On the other hand for a/ ~ 0.98 few 4- 
cliques emerge while 5-cliques appear around ctf ~ 0.73. This means that the social 
networks is mainly composed by 3-cliques, i.e. agents sharing high mutual affinities, 
that are connected together to form larger cliques, for instance 4 and 5-cliques by 
weaker links, i.e. whose mutual affinities are lower than the above ones. 

Results reported in left panel of Fig. [5] show that for specific parameter values, 
still falling into the class deputed to the consensus dynamics, the model does not 
present the weak ties phenomenon: 3, 4 and 5-cliques are all present at the same 
time for large values of ctf. This is an important point that will deserve future 
investigations. Let us observe here that the observed differences stem from the 
social temperature. 

3. Conclusion 

Social system and opinion dynamics models are intensively investigated within sim- 
plified mathematical schemes. One of such model is here revisited and analyzed. 
The evolution of the underlying network of connections, here emblematized by the 
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Fig. 5. Number of 3, 4 and 5-cliqucs in the social network once consensus has been achieved. 
Parameters are TV = 100, AO c = 0.5, ct c = 0.5. Right panel, a = 0.5, the network exhibits the 
weak ties property. Left panel, a = 0.1, the network does not display the weak ties phenomenon. 

mutual affinity score, is in particular studied. This is a dynamical quantity which 
adjusts all along the system evolution, as follows a complex coupling with the opin- 
ion variables. In other words, the embedding social structure is adaptively created 
and not a priori assigned, as it is customarily done. Starting from this setting, the 
model is solved analytically, under specific approximations. The functional depen- 
dence on time of the networks mean characteristics are consequently elucidated. The 
obtained solutions correlate with direct simulations, returning a satisfying agree- 
ment. Moreover, the structure of the social network is numerically monitored, via 
a set of classical indicators. Small world effects, as well weak ties connections, are 
found as an emerging property of the model. It is remarkable that such proper- 
ties, ubiquitous in nature, are spontaneously generated within a simple scenario 
which accounts for a minimal number of ingredients, in the context of a genuine 
self-consistent formulation. 

Appendix A. On the momenta evolution 

The aim of this section if to present and sketch the proof of the result used to study 
the evolution of the momenta of the affinity distribution. We refer the interested 
reader to [6| where a more detailed analysis is presented in a general setting. 

For the sake of simplicity, let us label the N(N — 1) affinities values by 
xi, upon assigning a specific re-ordering of the entries. Hence x is a vector with 
M = N(N — 1) elements. As previously recalled ([5]), we assume each xi to obey 
a first order differential equation of the logistic type, once time has been rescaled, 
namely: 
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The initial conditions will be denoted as x®. 



Let us observe that each component xi evolves independently from the other. 
We can hence imagine to deal with M replicas of a process ruled by by (jA.lj) whose 
initial conditions are distributed according to some given function. We are interested 
in computing the momenta of the x variable as functions of the initial distribution. 
The TO-th momentum is given by: 

< > it) = (xiW)m+ M +( " M(t)r , (A.2) 

and its time evolution is straightforwardly obtained deriving (|A.2|) and making use 
of Eq. (TO) : 



d <r T m -> (t) - 1 d%T_ _ m ^ m -idxi 
dt < M 2-/ dt ~ 1 dt 

i=l 1=1 



N 

= Jfi E x ?~ lx i( l ~xi)=m (< x m > - < x m+1 >) . (A.3) 
i=i 

To solve this equation we introduce the time dependent moment generating junc- 
tion, G(£, t), 

oo 

G(U)-= ^r<x m >(t). (A.4) 

m— 1 

This is a formal power series whose Taylor coefficients are the momenta of the 
distribution that we are willing to reconstruct, task that can be accomplished using 
the following relation: 

1 r) m C 

<x m >(«):= . (A.5) 

m! a£ m {=o 

By exploiting the evolution's law for each xi, we shall here obtain a partial 
differential equation governing the behavior of G. Knowing G will eventually enables 
us to calculate any sought momentum via multiple differentiation with respect to £ 
as stated in (|A.5|) . 

On the other hand, by differentiating (|A.4j) with respect to time, one obtains : 

^ = E ^ dt = E m ^ m (< xTn > - < xm+1 >) > ( A - 6 ) 

m > 1 m > 1 

where used has been made of Eq. (|A.3|) . We can now re-order the terms so to 
express the right hand side as a function of G El and finally obtain the following 
non-homogeneous linear partial differential equation: 

d t G - (£ - l)d £ G = | . (A.7) 

c Here the following algebraic relations are being used: 

£0 £ G(£,t) = £ m ? m_1 < x m >= m Z m < xm >' 
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Such an equation can be solved for £ close to zero (as in the end of the procedure 
we shall be interested in evaluating the derivatives at £ = 0, see Eq. (|A.5[) ) and for 
all positive t. To this end we shall specify the initial datum: 

G(£,0)= ££ m <* m >(0) =$(£), (A.8) 

m>l 

i.e. the initial momenta or their distribution. 

Before turning to solve (|A.7[) ■ we first simplify it by introducing 

G = e 9 namely g — log G , (A.9) 

then for any derivative we have 9»G = Gd*g, where * = £ or * = t, thus (|A.7[) is 
equivalent to 

dtg - (f - l)% = ~ , (A.10) 

with the initial datum 

5 (£,0) = ^)=log$(£). (A.ll) 

This latter equation can be solved using the method of the characteristics, here 
represented by: 

! = -(*-!), (A.12) 
which are explicitly integrated to give: 

£(i) = l + (£(0)-l)e- 4 , (A.13) 
where £(0) denotes at t = 0. Then the function u(£(t),t) defined by: 

mtU):= m 0)) + J* i + m ]_ i)e _ s d s, (A.14) 

is the solution of (|A.10[) . restricted to the characteristics. Observe that u(£(0),0) = 
</>(£(0)), so (|A.14j) solves also the initial value problem. 

Finally the solution of (jA.llj) is obtained from u by reversing the relation be- 
tween and £(0), i.e. £(0) = (f (i) - l)e* + 1: 

g(£,t) = cj>((Z-l)e t + l)+\(Z,t), (A.15) 

and 

^^iMl = ^ j2 c m_1 < * m >= ? ( m - 1 )^ m ~ 2 < xm > 

^ m>l m>l 

= J2 (™- l)?" 1 " 1 < z m > 

m>l 

Renaming the summation index, m — 1 — > m, one finally gets (note the sum still begins with 
m = 1): 

=J2rr>4 m < x m+1 > . 



I 



I 



February 3, 2010 10:5 WSPC/INSTRUCTION FILE CarlettiFanelliRighi 



14 T. Carletti, D. Fanelli and S. Right 



where A(£, t) is the value of the integral in the right hand side of (| A. 141) . 

This integral can be straightforwardly computed as follows (use the change of 
variable z = e _s ): 



A = 



which implies 



A = - 



1 



l + (£(0)-l)e- 



ds = 



—dz 



1 



dz 



1 



m - 1 



l x z l + (£(0)-l)* 

i + logCl + ^Oj-lje-O-log^O) 



z l + (£(0)-l)«' 
= -logz + log(l + (£(0)-l)z) 



(A.16) 



According to (|A.15j) the solution g is then 

t) = <j>((£-l)e t + l)+t + log £ - log((C - l)e 4 + 1) 
from which G straightforwardly follows: 

G(£,t) = $ ((£ - l)e 4 + 1) ^' 



(A.17) 
(A.18) 

(A.19) 



(£-l)e* + l' 

As anticipated, the function G makes it possible to estimate any momentum 
(| A.5|) . As an example, the mean value correspond to setting m = 1, reads: 



< x > (t) = d s G 



+ d> (1 + (£ - l)e*) 



(1 + (£ - l)e 4 ) e* 



(£-l)e* + l 
(e-iy + l-^e 4 - 



(l + (e-l)e*r 



1-e 



T $(l-e 4 ). 



(A.20) 



In the following section we shall turn to considering a specific application in the 
case of uniformly distributed values of affinities. 



A.l. Uniform distributed initial conditions 

The initial data x® are assumed to span uniformly the bound interval [0,0.5], thus 
the probability distribution ip(x) clearly reads o 



ip(x) 



2 if x G [0, 1/2] 
otherwise 



and consequently the initial momenta are: 

nl r l/2 



< x m > (0) 



o Jo 



1 1 

m+ 1 2™ 



(A.21) 



(A.22) 



d We hereby assume to sample over a large collection of independent replica of the system under 
scrutiny (M is large). Under this hypothesis one can safely adopt a continuous approximation for 
the distribution of allowed initial data. Conversely, if the number of realizations is small, finite 
size corrections need to be included [6]. 
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On the evolution of a social network 15 

Hence the function $ as denned in (|A.8|) takes the form: 

$(£)=V S_. (A.23) 

KS> ^ m+ 1 2 m v ' 

m>l 

A straightforward algebraic manipulation allows us to re-write (|A.23[) as follows: 



y m 1 TV- mj 1 r z , 1 1 



_ , / J2 zmdz = ~ 1 d* = -l--log(l-y), (A.24) 

yj ^ 2/J 1-z y 



*(0 = -l-|log(l-| ] . (A.25) 



thus 



We can now compute the time dependent moment generating function, G(£,i), 
given by (|A.19|) as: 



(C - l)e* + 1 
and thus recalling (|A.5[) we get 



(£-l)e* + l 

'-(^TFTT 10 ^ 1 



(A.26) 



<x> v = — ~V^r l0g {—) (A - 27) 

2 e 2t 4e 2t , /e* + lA 2e 2 



( e t_l)2 ( e t_i)3 6^ 2 J (e*-l) 2 (e* + l) ' 

Let us observe that < x > (t) deviates from the logistic growth to which all the 
single variable Xi(t) does obey. 

For large enough times, the distribution of the variable outputs is in fact con- 
centrated around the asymptotic value 1 with an associated variance (calculated 
from the above momenta) which decreases monotonously with time. 

Let us observe that a naive approach would suggest interpolating the averaged 
numerical profile with a solution of the logistic model whose initial datum x° acts 
as a free parameter to be adjusted to its best fitted value: as it is proven in [6] this 
procedure yields a significant discrepancy, which could be possibly misinterpreted 
as a failure of the underlying logistic evolution law. For this reason, and to avoid 
drawing erroneous conclusions when ensemble averages are computed, attention has 
to be payed on the role of initial conditions. 
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